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We show that a standard calculation of the splitting func- 
tion for the nonsinglet structure function hi does not lead 
to the expected result Q. The calculation is compared to 
the corresponding derivation of the splitting function for the 
nonsinglet polarized structure function g\. We analyze pos- 
sible explanations for the unexpected result and discuss its 
implications. 

PACS: ll.10.Hi; 11.40.-q; 11.55.Ds; 13.88.+C 



I. INTRODUCTION 

In the Drell-Yan process a quark- antiquark pair is an- 
nihilated to a virtual photon, which produces lcpton 
pairs. The corresponding cross section contains a contri- 
bution from quarks of different chirality: a quark origi- 
nating from one nucleon may return to this nucleon in the 
squared scattering amplitude with its chirality flipped. 
This is possible because the intermediate photon does not 
remember the distribution of quark chirality on the scat- 
tering nucleon. This gives rise to the leading twist chiral 
odd distribution function hi , which first appeared in the 
discussion of the transverse polarized Drell-Yan process 
0-pl- One should emphasize that such chiral odd quark 
distribution functions do not appear in totally inclusive 
deep inelastic scattering because the quark chirality is 
conserved at each vertex of the scattering process. 

The quark helicity distribution g\ has a partonic in- 
terpretation in the chiral basis. This is not the case for 
h\ which gets a partonic interpretation only by switch- 
ing to the transversity basis j(| . Looking on transversely 
polarized nucleons h\ is interpreted as the probabilistic 
asymmetry of finding quarks in eigenstates of the trans- 
verse Pauli-Lubanski vector with eigenvalue ±1/2. One 
should be aware, that g\ looses its partonic interpretation 
in the transversity basis. 

There are several attempts to measure the transver- 
sity distribution h%, see e.g. For general reasons 
the transversity distribution has to obey restrictions as 
Soffers inequality for twist-2 distribution functions 
2\hi(x)\ <h{x)+gi(x). _ _ 

The anomalous dimension of hi and the correspond- 
ing splitting function can be found in literature |l|,|l2|. 
This result is found by direct evaluation of the splitting 
function. Nevertheless the behavior at Bjorken-a; equal 
to 1 remains unfixed in this procedure. In the case of the 
helicity distribution gi the behavior for x = 1 is fixed by 
a corresponding sum rule. For hi such a sum rule is lack- 



ing. So the full hi splitting function was determined by a 
general argument, which demands the splitting function 
of hi to behave exactly as the splitting function of gi for 
x | 1. To our knowledge, no explicit calculation (without 
additional arguments) of the full splitting function of the 
transversity distribution exists until now. 

This paper is organized as follows. First, the forward 
scattering amplitudes corresponding to the helicity and 
transversity operator are defined. We derive the leading 
order anomalous dimension of gi and hi using a disper- 
sion relation and show that the behavior of the splitting 
function for hi does not match the expected result found 
in literature in the limit x | 1. We discuss possible ex- 
planations and implications of our result. 



II. DEFINITIONS 



A. Helicity operator 



The forward scattering amplitude 



(p, s\ r \j„ 5 (y)ju(o)+jAy)j^m \P, s) (1) 



contributes to the nonsinglet polarized structure func- 
tion gi(x,Q 2 ) := gl{x,Q 2 ) - g[ l (x,Q 2 ), where 
3v{y) = ^iy)l^{y) is a vector current and j^{y) = 
ip(y)jfij5ip(y) denotes an axial- vector current. gi{x,Q 2 ) 



depends on all possible Lorentz invariants Q 



and 



Bjorken x with x — Q 2 /2p ■ q. Here p and q are the four- 
momenta of the hadron and virtual photon, respectively. 
The latter defines the momentum scale of the scattering 
process. T denotes the time-ordered product and sum- 
mation over flavor indices is assumed. In order to get 
the part relevant for gi one has to project the forward 
scattering amplitude on the correct Lorentz structure: 



,9i 



1£ livXpQ 

p ■ q 



l 



Abs Tjjj, 



(2) 



v ■ q 



Abs denotes the absorptive part, i.e. the discontinuity 
accross a cut fl3|| . 

If higher twist operators are neglected, the moments of 
the polarized structure function can be identified with the 
twist-2 operator matrix elements A giJ i in the framework 
of the operator product expansion: 



1 



M gun = I dx x n g 1 {x) 



(3) 



Note, that we dropped the Q 2 - and the renormalization 
scale dependence, which are not important for our pur- 
pose. C n is the Wilson-coefficient which depends on the 
renormalization scale and the strong coupling. It con- 
tains the full perturbation series of leading twist. For ex- 
ample the coefficient for n — corresponds to the Bjorken 
sum rule. 



B. Transversity operator 

Let us define the following forward scattering ampli- 
tude (i 



14-lJ): 



T ll {p,q,s) = -j d 4 ye^x 

(p,s\T Mt/)js(0) + js(y)j^(0)} \p, s) ■ (4) 

jfj,5{y) is a nonsinglet axial- vector current as for g\ and 
js(y) = ' l P(y)' l P(y) a nonsinglet scalar current. Again, 
summation over flavor indices is assumed. This opera- 
tor defines a structure function, which describes vertices 
without conservation of chirality. Such vertices do not ex- 
ist in QCD. Nevertheless, quark chirality may be flipped 
in Drell-Yan processes, in which quarks of different chiral- 
ity may be connected to one virtual photon. This process 
is not suppressed by powers of Q 2 for large momentum 
transfers such that this contribution to the cross section 
is of leading twist. Consequently, the structure function 
hi defined by Eq. (Q) may be connected to the tranver- 
sity distribution. Indeed, it was shown on leading twist 
level that hi can be identified with the transversity dis- 
tribution (cf. the Appendix). 

As for the polarized structure function, we get the 
contributions of the forward scattering amplitude to the 
structure function hi by projection on the corresponding 
Lorentz structure: 



hi(x)s±^ = - AbsT) 



The corresponding moments are defined as in || 
l 

Mx,» = J dxx n [hi(x) + (-) n hi(-x)] 
o 



(5) 



(0) 



The last term again corresponds to the leading twist con- 
tribution in the operator product expansion with the Wil- 
son coefficient C±_ n and the chiral odd operator matrix 
element A\ „. 



III. DERIVATION OF SPLITTING FUNCTIONS 

A. Wick expansion 

The splitting functions for the operators in Eqs. (jl]) 
and (^) are determined by the divergent part of the for- 
ward scattering amplitudes. By Wick expansion of the 
current products we find 32 Feynman graphs contributing 
to the splitting functions. These diagrams can be divided 
into two classes which transform into each other by sub- 
stituting the gluon momentum k — > —k. Since the gluon 
loop integrals are symmetric it is sufficient to consider 16 
graphs. The forward scattering amplitude T is real on 
the real axis, i.e. T* = T . The remaining 16 graphs split 
up into 8 diagrams and their complex conjugates - which 
contain the same loop integral and a Lorentz structure 
in inverse order. The 8 relevant diagrams left consist 
of 4 diagrams (given below) with different loop integrals 
and the corresponding exchange graphs. The latter are 
obtained by substituting the momentum of the virtual 
photon q — » — q and interchanging the external vertices 
simultaneously. The Feynman diagrams are depicted in 
Fig.pTA] 
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FIG. 1. The Feynman graphs contributing to the forward 
scattering amplitude with T u = 1,7,,, respectively: The upper 
graphs represent the selfenergy (left) and the boxgraph (right) 
and the lower ones depict the vertex corrections. 

The amplitude associated to the selfenergy reads 

1 



M s = -C F u(p)T v 
1 



— Iplpl x 
irj 



i>+ i + irj 



— Iy.l5u{p) I] 





S ' 



(7) 



where Ig is given in the Appendix. In analogy we get for 
the vertex corrections 

_ 1 

Mvi = C F u{p)^p-f a Tvll3l P -T— 

p + 

and 



IT) 



(8) 



2 



M V 2 = -C F u(p)T v - 



17] 

a/3 



(9) 



and finally for the boxgraph 



M B = C F u(p)jp-f a T I/ -f l 3j tl -f 5 -f\j p u(p) Ig PX . (10) 

The integrals I^f and I c jf X can be found in the Appendix 
too. The exchange graphs are not shown. r„ stands for 
7„ or for identity accounting for the helicity or transver- 
sity operator, respectively. We have to isolate the diver- 
gent parts contributing to the Lorentz structures which 
correspond to g\ or hi. 



B. Dispersion relation 

In order to calculate the contributions of the above 
Feynman diagrams to the splitting functions we have to 
consider terms like 



1 



17] 



p'+d+ir] Q 2 {uj -I) +rf 



(11) 



with 



(p + q) 2 = -Q 2 (1-lo) and w = l/x. (12) 



The limit t] j cannot be taken easily at this stage of 
calculation. For example the Dirac identity 



lim_!_ = p(I ) Ti7T S( X ) 



(13) 



T denotes the forward scattering amplitude accounting 
for <?! and hi, i.e. T Miy and T^, respectively. These ampli- 
tudes are analytic functions of Q 2 and x. They possess 
a pole at x = 1 and a cut on the real axis for x < 1. 
Through the absorptive part Abs T we isolate the discon- 
tinuity accross the cut. The dispersion relation Eq. (|l5| ) 
connects the behavior of T in the regular (but unphys- 
ical) region x > 1 with the physical region x < 1, i.e. 
the coefficients of the Laurent expansion for large x cor- 
respond to the moments of the structure functions. Hav- 
ing expanded the amplitude T, the full information about 
the singularity structure is contained in the correspond- 
ing coefficients. In this representation the limit 7] J, 
can be taken for each coefficient separately. In this way 
the anomalous dimensions are calculated. By taking the 
inverse Mellin transformation we find the corresponding 
splitting functions. 

In the case of gi crossing symmetry holds, so that 
Eq. (|l5| ) gives rise to even moments only. Odd moments 
cancel identically and there is no way of getting informa- 
tion on the missing moments. Nevertheless, in order to 
calculate the splitting function we will have to evaluate 
the inverse Mellin transformation of the anomalous di- 
mension, which, for this reason should be known for all 
moments. As it will become evident later, our main in- 
terest lies on the singularity structure for x f 1, in partic- 
ular we are interested in contributions like £(1 — x) to the 
forward scattering amplitude T. The (5-distribution is el- 
ement of the Schwartz-Sobolev space S_i ([0, 1]), cf. Jl7| . 
Its domain of definition is the dual space &i ([0, 1]), i.e. 
{S x \f) = f{x) exists for all functions / that are element 
of 6 1 ([0, 1]) C L 2 ([0, 1]). This is also the necessary re- 
quirement for the definition of the plus-function 



is not sufficient to deduce the resulting distributions. Al- 
though the imaginary part of the r.h.s. of Eq. ( |Tl| ) gen- 
erates some (5-type distribution 



((i-x); 1 |/)= / 





/(*) - /(I) 



(16) 



hm — — — 

nio Q 2 (u - 1) + 7] 



-iTtd 



OJ-l 



= -2in^Q 2 (u-l) 8{Q 2 {w - 1)) 

(14) 

this distribution does not contribute to the moments in 
Eqs. (§) and (||) owing to the measure dx. 

We deduce the singularity structure of the splitting 
functions using the well-known dispersion relation (see 
e.g. PI) 



00 f 1 

W J 



T(y) 



ta 



n+l 



dy y n x 



which will be of importance later. The singular behaviour 
of the splitting function at x | 1 corresponds to high 
moments. But for large n the weight functions x n form a 
Cauchy series with respect to the norm of the Schwartz- 
Sobolev space 6i ([0, 1]) 



x n \\k=0 



(17) 



The convergence is even stronger (0(l/n 3 )) for the L 2 - 
norm, but, e.g., it does not persist for the 62-norm. This 
Schwartz-Sobolev space 6 2 ([0, 1]) is the domain of def- 
inition of the derivative of the ^-distribution (S'\f) = 
—f'{x)- If we now assume that the singularity structure 
of the splitting function does not contain such distribu- 
tions, but only the Dirac-5 or the plus-function, we may 
restrict ourselves to the (lowest) Schwartz-Sobolev space 
61 ([0, 1]). Then the inequality of Schwarz 



- (AbsT(y) + (-l) n AbsT(-y)) . (15) 

7T 



<<?!/} < Halle.! 



(18) 



3 



implies that 



dxx n+1 AbsT 



dxx n AbsT -> 



(19) 



holds for n 1 oo. Note that the contribution of 5(1 — x) to 
the moments is a constant. From Eq. ( |l9| ) the behavior 
of the anomalous dimension for n f oo can be determined 
by the even moments only. The singular part T S j ng of T 
corresponds to high moments, so that it is already fixed 
in this way (as long as it can be defined within 6_i). 

Now let us analyze the regular part T reg of the forward 
scattering amplitude T = T reg +T s i ng . Since T is assumed 
to be an analytic function, we suppose that its moments 
possess an analytic continuation into the complex plane, 
i.e. we presume that 



M(s) 



dxx s AbsT r , 



(20) 



We have computed all one-gluon-exchange diagrams 
mentioned above with the vertex T v =~f v . In the frame- 
work of dimensional regularization with d = 4 — e the 
divergent part appears as coefficient of 2/e. In order 
to analyze the calculation we present the results for all 
graphs separately. After projection onto the appropriate 
Lorentz structure 



1 //// 



2ie llu x p q A s 
P ■ Q 



we get using the results from the Appendix 



2 ^ 



,71+1 



(23) 



(24) 



Note that this geometrical series can be summed up to 
yield — ui) in the unphysical regime ui < 1. For 

the physically interesting absorptive part a 5(1 — x)- 
distribution is generated. The results for the vertex cor- 
rections and the boxgraph read 



is an analytic function for s > 1. This assumption is very 
reasonable but not necessary for our main result, which 
mainly concerns the singular part. After a conformal 
mapping of the complex plane onto the unit sphere §2 
the known even moments describe a convergent series of 
the (supposed to be) analytic function M(s). Since we 
have subtracted the singular part - which is completely 
determined by the behaviour of the high moments - the 
remaining regular part approaches zero in this limit 



dx x n Abs T re 



(21) 



Now the theorem of identity for analytic functions can 
be applied. (If two analytic functions coincide at a con- 
vergent series, then they coincide everywhere.) It follows 
that the result for even moments remains correct for odd 
moments as well. This means that we get all moments by 
calculating the even moments only (under the assump- 
tions made). 

We require the same assumptions to hold for h\ and 
proceed in the same way. 



2 

£ n=0 



n+1 i+E 



k=l 



k + 1 



, (25) 



M 



9 i 9 \ 

° n=0 V fe=l * / 



2 2 

Mb,h* = -Cfv-E^TH — ^ — \ 
e (1 + n)(2 + n) 

n— ' 



(27) 



The gluon propagator was taken in Feynman-gauge so 
that it is reduced to the g M1/ -term. Nevertheless, we 
checked that the sum of the contributions arising from 
the fc^fc^-term vanishes as required by gauge invariance. 
Taking these results together we find 



4 2 



1 + n (l + n)(2 + n) 



+ 45 n - 3 , (28) 



C. Results for gi 

The forward Compton scattering amplitude corre- 
sponding to the polarized structure function gi reads af- 
ter expansion in u> = 1/x: 



T 



A — ^^fiuXp 9l - n 9l , n 



v ■ q 



(22) 



where 



n 1 



k=l 



(29) 



The leading order anomalous dimension becomes 
C F ( 4 2 



7gi,n 



(4tt) 2 \1 + n (l+n)(2 + n) 



45„ - 3 , 



(30) 



and the inverse Mellin transformation generates 



4 



P g > (*) = $p (4 - 2(1 3ff(l-*)) 

(31) 

for the splitting function. This result is well known in 
literature Q. Usually, the distribution 35(1 — x) is fixed 
by the sum rule 



dxP gi (x) = 



(32) 



and is not calculated directlyfas we did here by using 
the dispersion relation Eq. (jl5|)). Due to conservation 
of helicity the integral over the quark distribution does 
not change with Q 2 , i.e. 7 91 ,o = 0. As transversity is not 
conserved a corresponding sum rule for the transversity 
distribution is not known. 



and after an inverse Mellin transformation to the splitting 
function: 



P x (x) = 



C F 
(4tt) 5 



(1-x), 



(40) 



This is the splitting function corresponding to the op- 
erator defined in Eq. (|). Additionally, according to 
the method of Ioffe and Khodjamirian, this result may 
be identified with the splitting function of the (twist-2) 
transversity distribution as defined in Eq. ( [43| ) in the Ap- 
pendix. However, our result Eq. ( (40|) differs from the one 
found in literature for the transversity distribution 111] 



at x = 1. 



(4tt) 2 



(1-x). 



35(1 



(41) 



IV. DISCUSSION 



D. Results for hi 

In the same way as for the polarized structure function 
we expand Eq. (j4|) inw: 



n=0 



(33) 



In this case we find for the one gluon exchange graphs 
(again using the results from the Appendix) after projec- 
tion on s± fl 



(34) 



n=0 



2 / 2 \ 

n=0 V k=l / 



2 / 2 \ 

Mv2,» = C F s ± »- £ 1 + £ fe-T ' ( 36 ) 

n=0 \ fc=l / 



Mb,h = 0. 



(37) 



The boxgraph does not lead to any divergent contribution 
to the Lorentz structure sx^- 
Adding up all graphs we get 



if 4 



e * — ' \ 1 + n 

n—0 



AS, 



This leads to the anomalous dimension 
C F ( 4 



T-i,' 



(4-n-) 2 \1 +n 



4S r 



(38) 



(39) 



The splitting functions Eqs. ( |3l| ) and (40) for gi and 
hi, respectively, were calculated using the same method. 
In the case of the operator corresponding to the longitu- 
dinally polarized structure function the expected result 
was found. Nevertheless, the analogous calculation for 
the transversity operator leads to a result which differs 
from the one in literature. We want to elucidate possible 
reasons for this discrepancy. 



A. Sum rules 



As emphasized in Sec. [II B| a direct calculation of the 
split ting function without using the dispersion relation 
Eq. ( |15|) does not fix the behavior at x = 1. This is done 
by the sum rule Eq. ( |32| ) in the case of the helicity asym- 
metry operator. We showed above, that the dispersion 
relation leads to exactly the result expected from the sum 
rule without any further argument. 

One may observe that the splitting functions for the he- 
licity asymmetry and the transversity operator Eqs. ( |3l| ) 
and (^) differ by an additional term proportional to 1— a; 
which vanishes for x = 1 and an additional 5(1 —x). This 
distribution enters P gi (x) with the same weight as the 
missing term found in literature Eq. (41) for P±. This is 
not accidental. As for the transversity operator no sum 
rule corresponding to Eq. ( |3^ ) is known, the behavior for 
x = 1 of the transversity splitting function was fixed with 
the general argument ]j], |l9] , p0| , that the probability of a 
transversity flip should vanish for x — 1. It follows from 
this (reasonable) requirement that the splitting functions 
P± and P gi should behave identically at x = 1. We want 
to emphasize that by using this argument, the divergent 
structure at x = 1 is not calculated explicitly but fixed 
by hand. 



5 



B. Comparing both calculations 

The Feynman graphs contributing to the splitting 
functions differ solely by the vertex r„ which is 7^ in 
the case of the helicity asymmetry operator and unity 
for the transversity operator. For this reason the diver- 
gent part of one gluon exchange contributions are equal 
for the self energy A4s as well as for the vertex correction 
A4v2 which do not involve the differing vertex into the 
loop integral. Comparing the explicitly calculated results 
Eqs. (U) and (ffj) with Eqs. @ and ©, respectively, 
this statement is confirmed. 

Looking at the results for the boxgraph Eqs. ( p7|) and 
( |37| ) we observe a difference. The boxgraph is not diver- 
gent for the transversity operator due to an extra fac- 
tor 4 — d appearing by projection on the corresponding 
Lorentz structure. After an inverse Mellin transforma- 
tion the contribution Eq. ( p7| ) to the anomalous dimen- 
sion is translated into the 1 — x term in Eq. (|3l]). This 
term vanishes for x = 1 and for this reason is not a good 
candidate to explain the difference at x = 1. 

In this way only the vertex correction M.vi with a 
gluon exchanged at the j v or the unity vertex, respec- 
tively, may be responsible for the difference at x = 1. 
Comparing Eqs. ( p5f ) and (35) the result differ by an 
additional constant contribution to the anomalous di- 
mension. An inverse Mellin transformation of a constant 
leads exactly to the additional term — 3<5 ( 1 — x) for the 
polarized splitting function P gi . It is plausible that two 
graphs with one gluon exchanged at different vertices do 
not lead to the same result. Anyhow, we will have to 
look for a conceptual problem if we want to save the gen- 
eral statement of vanishing transversity flip probability 
at x = I. 



C. Assumptions 

We based the identification of the calculated anoma- 
lous dimension with the one of the transversity distri- 
bution on the argument of Ioffe and Khodjamirian p[ . 
Here it was shown, that the operator used in the for- 
ward scattering amplitude Eq. ([I|) indeed has a relation 
to the definition of the transversity distribution Eq. (43). 
Their argument remains correct on the leading twist- 
level, which is appropriate to our purpose. The identifi- 
cation may be affected by anomalies. But these appear 
first in second order in the strong coupling. This implies 
that our first order calculation remains unaffected from 
possible anomalies. 

In addition the method of Ioffe and Khodjamirian is 
based on light cone dominance. One may contest, that 
this assumption is not appropriate to describe the split- 
ting function at x — 1. Here a potentially unbounded 
number of zero- momentum gluons appears, which are 
not well described by free quark operators. Neverthe- 
less, the identification of the structure function (defined 




)) with the transversity distribution (defined in 
ff) on the leading twist level and in leading or- 
der of perturbation theory is unaffected by the above 
consideration, so that the leading order anomalous di- 
mension should be reproduced correctly. Consequently, 
if the method of Ioffe and Khodjamirian is applicable to 
the problem under consideration, we will have to look for 
other reasons for the discrepancy. 

Our calculation of the anomalous dimension for g\ and 
h\ rely both on the assumption that Abs T in Eq. (|l5|) is 
element of the Schwartz-Sobolev space &-\. The result 
Eq. (|3^) for gi is in accordance with literature which is 
not the case for h%, so that the above assumption may be 
questioned for the transversity operator. If AbsT M is not 
element of the Schwartz-Sobolev space <3_i, the irregular 
part contains other contributions than a ^-distribution. 
For example derivatives of <5 could appear. As the dif- 
ference of the results Eqs. ( f4(i| ) and is element of 
6_i, it cannot be explained in this way. This means 
that the reason for the discrepancy is not related to our 
assumptions. 

The assumption that the regular part of the forward 
scattering amplitude is an analytic function, is com- 
pletely independent of the problem at x = 1 . Analyticity 
is necessary to ensure that our calculation of the terms 
with x < 1 reproduces the correct regular part of the 
splitting function. 



D. Implications 

It follows from the analysis of the assumptions that our 
result for the transversity splitting function can be inter- 
preted in two alternative ways: First, the correspondance 
of Eqs. (||) and ( p3| ) shown by Ioffe and Khodjamirian 
contains a conceptual difficulty at x = I. Alternatively, 
our result implies a non- vanishing probability of transver- 
ity flip by emitting infrared gluons from quarks. This is 
really hard to believe. Especially, the number of infrared 
gluon emissions is arbitrary large so that the transver- 
sity becomes completely statistically distributed. On the 
other hand the regularization of the infrared regime was 
done by hand, exactly in order to avoid this surprising 
interpretation. It was not deduced from first principles 
of QCD. Our calculation shows that (under the assump- 
tions we made) the vanishing transversity flip probability 
for zero-momentum gluons does not follow analytically. 



APPENDIX 

A. Distribution functions 

The polarized quark distribution is defined in leading 
twist by 



G 



— e lXx (p, s\ V'(0)7 P 7 5 V'(An) \p, s) = gi(x, q 2 )p^s ■ n . 

(42) 

Again x = —q 2 /(2p ■ q) is the Bjorken variable. n u is a 



light cone vector with n 2 
the proton four-momentum with p 2 



of dimension (mass) , p is 
m 2 and p ■ n = 1 , 



and s is the proton spin four-vector with s 



-1. In 



deep inelastic electron proton scattering the distribution 
gi describes the asymmetry of nuclcons polarized parallel 
or antiparallcl to the longitudinally polarized lepton. 

The definition of hi in terms of an operator matrix 
element reads f0Ol: 



is obviously infrared divergent only, which leads to a max- 
imum divergence of 1/e. In the case of the transversal 
structure function h± the boxgraph is supplied with a 
prefactor of 4 — d resulting from the according Lorentz 
structure. Therefore it yields no contributions to the 
anomalous dimension. On the other hand the boxgraph 
leads to a term ~ (1 — ar) (see Eq. (27)) for the polarized 
structure function g\. However, since we are mainly in- 
terested in the behavior of the splitting functions in the 
limit x 1 1, the boxgraph is not relevant for our purposes. 
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where the transverse part of the spin vector is defined by 
the decomposition s M = (s • n)p^ + (s • p)n^ + sx^. 



B. Integrals 

The loop integrals neccessary for the calculation of 
the Feynman diagrams in Sec. HI possess infrared as well 
as ultraviolet divergences, which have to be regularized. 
For our purposes we used the dimensional rcgularization 
scheme with d = 4 — e. Since we are interested in the 
anomalous dimension and therefore in the divergent parts 
only, we can evaluate the integrals in the e — > limit. For 
the integral appearing in the selfenergy graph we find 



d d k 



2,.4-d 



(p + q- ky- 



rrj 



(2-K) d [(p + q- k) 2 + rf] (k 2 + irj) 
= -- £ l(p + q) a -ir 1 }(l + 0(e)). 

Both vertex corrections involve the same integral 

d d k 



(44) 



^=<(47T)V 



2,,4-d 



(2tt) 



[(p — k) a — ii]][(p + q — k) 13 — irj] 
[(p + q-k) 2 + r; 2 ] \(p -k) 2 + rf] (k 2 + irj) 

= — g —(l + 0(s)). (45) 
e 4 

The integral corresponding to the boxgraph 



nr-w* j w 

[(p — k) a — irj\[(p + q — k) 13 — irj\[(p — k) x — irj\ 
[(p + q — k) 2 + r/ 2 ][(p — k) 2 + T] 2 ] 2 (k 2 + irj) 



2 ,,4-rf 



d d k 
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